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ABSTRACT 



We investigate the ideas that a sizable fraction of the interferometrically unresolved infrared emission of the nucleus of NGC 1068 
might originate from other processes than thermal dust emission from the torus. We examine the contribution of free-free or syn- 
chrotron emissions to the central mid- and near-IR parsec-scale emitting region of NGC 1068. Each mechanism is constrained with 
parsec scale radio data available for NGC 1068 in the 10 9 - 10" Hz regime, and compared to the highest-resolution interferometric 
data available in the mid-infrared. It is shown that the unresolved emission in the interferometric observation (SI pc) is still domi- 
nated by dust emission and not by contributions from synchrotron or free-free emission. As recent studies suggest, the interferometric 
observations prefer a clumpy structure of the dust distribution. Extrapolation of the radio free-free or synchrotron emission to the IR 
indicates that their contribution is <20% even for the unresolved fraction of the interferometric flux. 

The slope of the available radio data is consistent with a power law exponent a = 0.29 ± 0.07 which we interprete in terms of either 
free-free emission or synchrotron radiation from quasi-monochromatic electrons. We apply emission models for both mechanisms in 
order to obtain physical parameters. 

Furthermore, we aimed for quantifying the possible contribution of the accretion disk to the near-infrared emission. It has been 
suggested, that the unresolved A^-band flux in VLTI/VINCI interferometric observation at 46 m baseline (40% of the total A-band 
flux) might originate from the accretion disk. By using an accretion disk spectrum which was adjusted to the luminosity and black 
hole mass of NGC 1068, we found out that the expected accretion disk flux in the A"-band is negligible. Moreover, the scenario of 
detecting the accretion disk through holes in a clumpy torus is extremely unlikely. We conclude that all current IR data of NGC 1068 
trace the torus dust emission, favouring a clumpy torus. 

Key words. Galaxies: Seyfert - Galaxies: nuclei - Galaxies: individual: NGC 1068 - Infrared: galaxies - ISM: dust, extinction 



1. Introduction 

The nucleus of the nearby Seyfert 2 galaxy NGC 1068 (d ~ 
14.4 Mpc; a scale of 10 mas corresponds to ~0.7pc) is the prime 
target of the most recent high-spatial resolution studies of AGN 
in the wavelength r e gime from infrared (IR) to the radio (e.g. , 



Weigelt et al. 2004 yWittkowski et all 120041: UaffeetalJ l2004t 
Gallimor e et al.ll2004l) . For the first time, it is possible to build- 



up a genuine nuclear spectrum of the immediate vicinity of the 
AGN. Furthermore, near- (NIR) and mid-infrared (MIR) inter- 
ferometric observations resolved the nulcear IR emission so urce 
dWittkowski et al.lll998tlWeigelt et al.ll2004t iJaffe et al.ll2004l) . It 
is supposed that the nuclear radio emission comes from an accre- 
tion disk corona or an associated jet while the infrared emission 
is more likely caused by thermal dust emission in the circumnu- 
clear dusty torus (e.g.. lGallimore et al.|[J996l l2004t iHonig et all 
l2006h . 

It is, however, yet unclear how much each component con- 
tributes at different wavelength s . MI R VLTI/MIDI interferomet- 
ric observations by Jaffe et al. (2004) show that approximately 
20% of th e MIR emission originat es from an unresolved source 
;£ 0.7 pc. Wit tkowski et all d2004) report that K-band interfer- 
ometry with the VLTI/VINCI instrument at 46 m baseline re- 
veals a surprisingly high contribution of 40% to the whole nu- 
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clear flux by an unresolved source of ^ 0.3 pc. This can be ex 
plained by either radiation from a small compact source which is 
associated to the nuclear radio spectrum (e.g., Wittkowski et al. 
1998) or by substructure inside a clumpy dust torus (IHonig et al 
2006). 

In this paper, we use the combination of high-spatial reso- 
lution photometry, radiative-transfer simulations in the IR, and 
emission models in the radio to compare the extrapolated radio 
emission to measurements at IR wavelengths, in order to inves- 
tigate the origin of the resolved and unresolved the near- and 
mif-infrared emission. In Sect. [2] we summarize the observa- 
tional data used for our study. In Sect. [3] we describe the models 
which were used for the radio and IR regimes, respectively. In 
the following Sect. [4] we discuss the model results and analyse 
the maximum contributions of the radio component to the ob- 
served NIR and MIR emission of diffraction-limited and inter- 
ferometric data. In Sect. [3] the contribution of the accretion disk 
to the IR spectrum is quantified and compared to the interfero- 
metric VLTI/VINCI observations. In Sect. [6] we summarize our 
results. 



2. Observational data 

The location of NGC 1068's nucleus in the radio is assumed to 
be at the center of a 0.8 pc size, disk-like structure known a s 
SI, resolved at 5 and 8.4 GHz bv lGallimore"etail (fl997l 12004). 
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Table 1. High spatial resolution observations of the nulear emis- 
sion of NGC 1068 from the radio to the infrared regime. 



Frequency 


Flux spatial scale 
[Jy] [mas] 


Reference 


2.3 x 10 14 


0.0084 78 


this work 


1.8 xlO 14 


0.022 130 


this work 


1.4 xlO 14 


0.098 78 


this work 


1.2 xlO 14 


0.32 ± 0.03 78 


this work 


7.9 x 10 13 


1.5 700 


[1] 


6.7 x 10 13 


2.5 162 


this work 


3.8 x 10 13 


4.23 200 


[2] 


3.0 xlO 13 


3.33 200 


[2] 


2.5 x 10 13 


7.1 200 


[2] 


1.6 xlO 13 


9.6 290x180 


[3] 


1.2 xlO 13 


9.4 200 


[2] 


2.3x10" 


0.022 ± 0.008 1000 x 600 


[4] 


1.15 x 10" 


0.023 ± 0.05 1400 


[4] 


4.3 x 10 10 


0.019 ± 0.003 75 


[5] 


2.25 x 10 10 


0.016 ± 0.006 75 


[5] 


8.4 x 10 9 


0.0054 ± 0.0005 11.0x3.7 


[6] 


5.0 x 10" 


0.0091 ±0.0008 16.6x11.2 


[6] 


5.0 x 10 9 


0.0051 14x3 


[7] 


1.7 x 10 9 


< 0.0026 14 x 12 


[7] 


1.4 x 10 9 


< 6 x 10" 5 16.0 x 7.6 


[6] 


[1] Marco & Brooks (2003); [2] Bock et al. (2000); 


T31 Tomono et al. (2001); 141 Krios et al. (2006); 




T51 Muxlow et al. (1996); [61 Gallimore et al. (2004); 


T71 Rov et al 


(1998) 





In the mid-IR, a gaussian decomposition of VLTI/MIDI inter- 
ferometry data sets an upper limit of 0.7 pc x <1 pc on the siz e 
of the central source at wavelengths < 10pm (Jaffe et al. 2 0041). 
At 2 / mi, the combination of speckle obs ervations dWeigelt et alj 
2004b and VLTI/VINCI interferometry dWittkowski et al.l l2004) 
confirm the presence of a compact core with FWHM ~0.4pc. 
Shortward of 1 /urn, the NGC 1068 nucleus is fully obscured and 
undetected. 

The radio to IR to optical SED of N1068's nucleus is ex- 
tracted from subarcsec photometry currently available for this 
source. Overall, we sample spatial scales at radio, MIR and 
NIR wavelengths down to few parsecs. Radio data for the 
source SI, which is assumed to be NGC 1068 nucleus, is 
taken from the VLA/MER LIN 65 mas beam (Muxlow et al. 
1996; Gallimo re et al. 12004 ) and VLBA-observation dRov et al] 
1998; Gallimo reeraTll2004 . lGallimore et alj d2004l) report that 
they resolved the nuclear source SI at 5 and 8.4 GHz. SI ap- 
pears disk-like with a reference size of ~llmas (for detailed 
sizes see Tabl e [TJ). The milli meter data at 1 mm and 3 mm are 
taken from lKrips et alj (|2006) and correspond to a beam size of 
1 x 0.6arcsec. Because of the possible contribution of the jet in 
the 3 m m data, only the peak value as specified by Krip s et alj 
(2006) is considered. 

Data in the 8 — 25 pm range are taken from lBock et al.1 d2000t) 
and iTomono et alj d2001l) . using their reported fluxes in aper- 
ture diameters of ~200mas. This high resolution was obtained 
by deconvolving their 10 m Keck and 8.2 m Subaru images, re- 
spectively. It is important to note that both sets of deconvolved- 
extracted fluxes differ by <30%. 

Interferometry d ata from VLTI/MI DI in the 8-13 pm win- 
dow are taken from iJaffe et alj d2004|) . For our study, we used 
the correlated fluxes from the 78 m baseline which provides the 
highest resolution data currently available. The visibilities show 
that the silicate absorption feature becomes more pronounced 



with growing baseline. This was interpreted as an indication 
for clumpy substructure within the dust torus dJaffe et al] |2004; 
Honig et al. 2006). 5-20% of the total flux are unresolved at the 
78 m baseline. 

In the 1 - 5 pm rangqj we use our adaptive optics im- 
ages acquired with VLT/NACO, except for the L-band pho- 
tometry which was taken from Marco & Brooks (2003). The 
spatial resolution achieved in these observations varies from 
FWHM-78 mas in the 1.3 -2.44 pm range to FWHM-160 mas 
in the M-band. Nuclear fluxes were integrated in an aperture size 
comparable to those resolutions. These aperture sizes were se- 
lected as a compromise to minimize the contribution from the 
prominent IR diffuse emission surrounding the nucleus of NGC 
1068. The 1.3 pm measurement should be considered as an up- 
per limit. The corresponding frequency is the highest at which 
the core of NGC 1068 can be reliably distinguished from the 
strong extended emission. 

NGC 1068 is Compton thick in the X-ray regime. It shows 

S irominent extende d soft X-ray emission as seen by Chandra 
Young et alj|200ll) . The extended emission is preferentially dis- 
tributed along the ionization cone. XMM spectra show that most 
of this emission is resolved i n multiple emission line s and not 
much continuum is detected (Kin khabwala et al]|2002l) . Due to 
the Compton thickness of the source and extension of the emis- 
sion, we do not include X-ray photometry in our modelling of 
the nucleus. 



3. Models 

We focused our modelling efforts on two components: (1) the 
near- and mid-infrared component in the range of 10 13 - 10 14 Hz, 
and (2) the radio component at frequencies from 10 9 Hz to 
10 11 Hz. Our aim is to quantify the contribution of each com- 
ponent in different wavelength regimes. 

3.1. IR torus model 

For the infrared regim e of the NGC 106 8 SED, we used the 
clumpy torus models of lHonig et a l. (2006). With these models, 
we aimed at fitting the photometric data from 1 .2 - 23 x 10 13 Hz, 
simultaneously with the 78 m baseline data as observed with 
VLTI/MIDI. A grid of models has been calculated where we 
tested (a) different radial cloud distributions (radial distribution 
function oc r~ a with a - 1.0, 1.5,2.0), (b) inclincation angles, 
and (c) average numbers of (optically thick) clouds, No, in an 
equatorial line of sight. For each set of model parameters we 
simulated 10 different random arrangements of clouds. This ac- 
counts for the variations in SEDs which occur due to the ac- 
tual cloud distribution. Our best-fitting model is shown as blue- 
dashed lines in Figs.[TJ&[2] 

The modeled torus has a line-of-sight inclination of 90 deg, 
which corresponds to an edge-on view of the torus. In radial di- 
rection, the clouds are distributed according to a power-law of 
r~ 3 ^ 2 , which is well consistent with clumpy torus accretion sce- 
narios where the torus is located i n the sphere of influence o f the 
central supermassive black hole (Beckert & Duschlll2004l) . We 
assumed linear flaring throughout the torus, i.e. the height of the 
torus H increases with radius, H oc r. The half opening angle 



1 The near-IR wavebands are referred to as ./-band (1.3 pm <=> 2.3 x 
10 14 Hz), //-band (1.65 pm O 1.8 x 10 14 Hz), /f-band (2.2pm O 
1.4xl0 14 Hz), L-band (3.8 pm «■ 7.9xl0 13 Hz), and M-band (4.5pm «■ 
6.7 x 10 13 Hz). Mid-IR wavebands are referred to as iV-band (8 - 
13pm »~3x 10 13 Hz) and g-band (16 - 25 pm & ~ 1.5 x 10 13 Hz). 



S. F. Honig et al.: High-Spatial Resolution SED of NGC 1068 from Near-IR to Radio 



icr 



10 1 - 



10 u 



icr 1 - 



Synchrotron model 



* .«..—*"!" 



\o~" - / 

I 

/ 
/ 

lO" 4 - 



10 s 



1/ 






10' 




frequency i/ [Hz] 



Fig. 1. High-spatial resolution SED of the nucleus of NGC 1068, in comparison to thermal dust emission from a torus in the IR and 
non-thermal synchrotron emission in the radio. The blue line shows the clumpy torus model that was used in this study to model 
the single-dish IR photometry (diamonds) and interferometric observations by VLTI/MIDI (orange solid line; model: dashed-dotted 
blue line). The radio data (diamonds) are modelled by a spectrum F v oc y 1 / 3 exp(— (v/Vff) -21 ) (Synchrotron emission + free-free 
absorption; green-dotted line) and extrapolated into the infrared. The dashed line shows the synchrotron spectrum attenuated by 
dust obscuration (Ty = 16.8) to reproduce the /-band data, regimes. 



Table 2. Parameters of the synchrotron models A, B, and C for the nuclear source SI of NGC 1068: The linear scale s is the half 
thickness in model A & B and the radius of the emission spheres in model C. For the mean Lorentz factor (y) of the electron 
distributions (width Ay/ (y) = 2) there are only lower limits derived form the assumed upper cut-off frequency v c — 3 x 10 14 Hz of 
the synchrontron spectrum. The values for B M ; n is the minimum energy B-field for the claimed size of s - 5.6 mas in model A and 
upper limits for B and C. The electron density in the last column corresponds to Bjviin- 



Model s 

[mas] 



A 
B 



5.6 

> 2 x 1(T 7 

> 1.6 x 1Q- 3 



sxd 
[pc] 



<y> 



0.39 

> 1.3 x 10- 8 

> 1.1 x 10~ 4 



> 1.5 x 10 5 

> 1.3 x 10 4 

> 8 x 10 3 



[G] 



1.6 x 10~ 2 

<2.2 

<6.1 



n, 
[cm" 3 ] 



1.3 x 10-" 

<3.0 

<350 



of our geometrical configuration is ^40deg. Although a general 
power-la w H oc r x could have been chosen, it has already been 
tested in iHonig et all d2006l 120071) and no better fit to the IR 
data was found then. Furthermore, theoretical considerations and 
observational studies of clumpy tori seem to favour linear flar- 
ing g eometries ove r stronger r-dependence s of the torus height 
dHonig & B eckert 2 0071 iPolletta et al.l2008l) . The obscuration of 
the nucleus in the presented torus model is characterised by an 
average of No - 10 optically thick dust clouds in an equato- 
rial line of sight. For the edge-on geometry that we modelled, 
this corresponds to average number of clouds, N, which ob- 
scure the AGN along the line of sight towards the observer (i.e., 
N - No - 10, as used in Sect.©. 

As compared to other model parameters in our grid, it si- 
multaneously reproduces (1) the NIR slope, (2) the depth of the 



silicate feature, (3) the slope from A-band to Q-band, and (4) 
the correlated MIDI flux. Due to the ambiguity in actual cloud 
arrangement, there is som e degeneracy in the model parameter 
space (compare results in iHonig et al] 120061) . especially when 
fitting the SED alone. On the other hand, the presented dust dis- 
tribution parameters are the only one for which we found good 
agreement with the SED properties (1-3) and the interferomet- 
ric data (4). As a general guideline, a > 1.5 produces bluer NIR 
colors and le ss Q-band flux, and vice versa for a < 1.5 (see also 
discussion in lPolletta et al]|2008l) . 

The observed and modeled diameter of the dust torus is ~20- 
30 mas (~1.4-2pc), depending on wavelength and orientation. 
To fit the observed fluxes, a bolometric luminosity of the central 
AGN of 2 x 10 45 erg/s has been used. We note that this model lu- 
minosity is larger than inferred from the measured fluxes assum- 
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ing isotropic emission. The reasons for this are actual anisotropic 
emission of the torus (type 2 edge-on emission is smaller than 
type 1 face-on emission) as well as the covering factor of the 
torus which is < 1 , so that not all of the AGN bolometric lumi- 
nosity is absorbed by the torus. Total discrepancies can be up to 
one order of magnitude. 

3.2. Radio models 

In the radio regime at frequencies lower than 1 THz, we used 
two different emission models to reproduce the observed high- 
spatial resolution data. In order to study the possible influence of 
the radio emission sources to the observed IR emission, we will 
assume that the model spectra reach the near-IR regime without 
a high-frequency break. In a further step, the model spectra are 
attenuated by dust extinction so that the model fluxes are consis- 
tent with the observed near-IR fluxes. For that, we use different 
optical depths, Ty, so that the observed /-band flux at 1.2 fim is 
reproduced only by the radio component (i.e., 100% flux con- 
tribution). This procedure will give us a very strong upper limit 
to the possible contribution of the radio source to the Zf-band 
data where high-spatial resolution observations and interferom- 
etry is available. We will compare the used tv values with the 
torus model and the estimated broad-line region obscuration in 
the next section. 



Synchrotron radiation. The radio data of NGC 1068 shows an 
inverted spectrum (a > 0) between 5-230 GHz with an appar- 
ent turnover at ~ 10 9 Hz. From a least-square fit to actual ra- 
dio detections we obtain a slope of a = 0.29 + 0.07. This is 
consistent with a power law index a mono = 1/3, which can 
be produced by optically thin synch rotron radiation comin g 
from quasi-monoenergetic electrons (Beckett & Duschl| [l997). 
In such a scenario, the emission spectrum can be approximated 
by 



,1/3 



exp(-v/v c ) 



(1) 



with an exponential turnover, v c , at high frequencies. For the 
purpose of estimating the maximum possible NIR contribution 
of the radio source, we assume that v c <: 300 THz, so that the 
actual high-frequency cutoff can be neglected hereafter. At low 
frequencies, a turnover can occur, e.g. if the synchrotron source 
is located within dense plasma leading to thermal free-free ab- 
sorption Tff oc v which produces a sh arp turnover towards 
low frequencies (e.g jBeckert et allll996l) . As a result, the syn- 
chrotron spectrum is modified by a factor exp (— (v/vff )~ 2 ' ' ), with 
the low-frequency turnover at Vff, so that the flux can be param- 
eterised by 



5 v <xv 1 / 3 exp(-(y/v ff )- 21 ) , 



(2) 



neglecting the exp(-v/v c ) turnover at high frequencies. 

In Fig. Q] we compare the synchrotron model SED with 
observed radio fluxes. The low-frequency free-free absorption 
turnover occurs at vg - 3 GHz which is consistent with the upper 
limit on detections around 1.5 GHz. The spectrum was attentu- 
ated by standard ISM dust absorption of optical depth Tv = 16.8 
to reproduce the /-band observation. 

An alternative explanation for the low-frequency turnover is 
synchrotron self-absorption (SSA) within the source. This would 
result in a spectral turnover from the optically thin F Y oc y 1 / 3 to 
the optically thick case F v oc y 2 , which is less steep than in case 
of free-free absorption as discussed above. The optically thick 



part behaves like a thermal Rayleigh- Jeans branch in contrast 
to the standard power-law electron distribution which produces 
F v oc y 5 / 2 (e.g., Pacholczvk 1970). Given the detection of SI 
in NGC 1068 at 5 GHz and the upper limits at -1.4 GHz, the 
turnover should occur inbetween. The observed lower limit for 
the power law index in this frequency range (a ^3.5) favours the 
free-free absorption scenario over SSA. 

In a quantitative analysis, we inverted the low- 
frequency model approximation to synchrotron emission 
and self- absorption from quasi-monoenergetic electrons 
(Beckett & Duschl 1997). We suppose that the minimum energy 
assumption holds for our emission region to link B-field and 
electron density, n e . The emission region is asummed to be a el- 
lipsoid with apparent axes r, s (volume 4/3nr 2 sd 3 , and distance 
d to NGC 1068; disk-like if s < r). In case the low-frequency 
turnover is due to SSA and the peak of the spectrum (end of the 
v^-part) can be determined from observations, it is generally 
possible to fix the main model parameters, i.e. the fi-field 
strength, the electron density n e , and the Lorentz factor y. Here, 
we don't know the peak of the spectrum and have upper limits 
to the self-absorption frequency since we assume that the steep 
low-frequency turnover is due to free-free absorption rather 
than SSA. Thus, only limits of the physical parameters can be 
determined. 

For reference, we use the radio flux of the partly resolved 
core SI at vi = 5 GHz from (iGallimore et all 12004) of F Vl - 
9.1 + 0.8 mJy. The size of the elliptical source is 16.6 and 
11.2 mas for major and minor axis respectively. In Table|2]we use 
the half axis values corresponding to radii in the source model. 

Here, we present three synchrotron models for SI. All are 
based on the v 1/,3 -power law spectrum as described above. In 
the first case (model A), we assume that the disk is resolved 
in the radial and vertical direction with a half major axis r - 
8.3 mas and a half minor axis s - 5.6 mas. The corresponding 
source parameters are given in Table [2] In model B, we assume 
that the actual disk thickness is smaller than the minor axis, s < 
5.6 mas, so that the observed vertical extension might be due to 
disk inclination or warping. In this case we derive a lower limit 
for s from the absence of SSA down to 3 GHz. Thus, a lower 
limit for the mean Lorentz factor y of the electron distribution is 
obtained, as well as an upper limit for the corresponding B-field 
and also for the density of relativistic electrons in the radiating 
plasma (values see Table [2]). 

The altenative model C is motivated by the orientation of the 
major axis of S 1 which is alm ost perpendicular t o the large scale 
radio jet in NGC 1068 dMuxlow et al.lll996t IGallimore et"ai] 
2004). SI may, therefore, be associated with a perpendicular 
shock front close to the base of an outflow or the corona of 
an accretion disk. The large size of SI and its low brightness 
temperature of Zsghz = 3.5 x 10 6 K argues against a large 
scale shock. Instead, a corona above an accretion disk may be 
heated in locallized regions through magnetic reconnection or 
small scale shock fronts. The corresponding emission regions 
are much smaller than the size of the underlying disk. In this 
model C, we assume N ~ 100 unresolved coronal emission 
regions which together form SI and give it the appearence of 
a thick disk. For simplicity we assume spherical emission re- 
gions and use a 3 times larger size than implied by the SSA 
limit to avoid mutual absorption by emission regions superim- 
posed on the sky. We obtain lower limits to the sphere radius s 
and to a mean Lorentz factor (y), and upper limits for B-field 
and relativistic electron density n e . The results are summarized 
in Table |2 In this picture, the relatively large size of SI as com- 
pared to the small emission regions might be due to electron scat- 
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Fig. 2. Same as Fig.Q] but with our free-free emission model in the radio. The free-free spectrum was attenuated by dust with Ty = 6 
to reproduce the /-band data. 



tering not directly connected to the small emission regions (e.g. 
surrounding medium) and provides only an upper size limit. 

For NGC 1068, the possibility of synchrotron emission 
from quasi-mono e nerge tic electrons was already explored by 
Witt kowski et al.l U998) but with less radio data available. 
Similar models have been suggested for other galactic nuclei like 
SgrA*dDuschl &Lesch 1994; Beckert & Duschl 1997b or_M81 
dReuter &Leschlll996l) . According to IWittkows ki et al.l (|2004|) . 
such an emission mechanism would require a source size of the 
order of 100AU and electron densities of 10 3 ctrT 3 . While the 
size is consistent with our estimates, we require lower electron 
densities in the presented analysis. 

Free-free emission. The fluxes at the highest radio frequencies 
> 100 GHz are consistent with a flat spectrum. This opens the 
possibility that the radio spectrum is produ ced by therm a l free- 
free e mission, as already suggested by Gallim ore et alj d 19971 
120041) . 

To investigate this possibility we setup a free-free emission 
model. The flux from thermal free-free emission at frequency v 
for a given electron temperature T is described by 



S Y = B v (v,T)(l-s- T «) 



(3) 



where B v (v,T) denotes the Planck function B v - Ihv^jc 1 ■ 
(exp(hv/(kT) - I)" 1 . The optical depth for the free-free process 
can be described as 



r ff = 1.13725 (l-e-S)ftKv,r) 

X (ZX 1/2 (JL)- 3 (JH.) 2 [L 

\K/ \GHz/ \cm 3 / \pc 



(4) 



Table 3. Parameters for the free-free emission model as de- 
scribed in the text. 



Model 


e/D 

[mas] 


I 
[pc] 


[cm" 3 ] 


T 

[K] 


Tv 


F 


11 


0.75 


8xl0 5 


1.3 x 10" 


6 



Here, the emission measure EM - J n^ dt is approximated by 
a slab with constant electron density n e and a typical size I. The 
function gg is the Gaunt factor which accounts for the differ- 
ence between semi-cl assical and q uantum mechanical treatment 
of free-free emission dKramerslll9 23; Gaunt 1930). For low fre- 
quencies where j ?ff > 1 we u se the thermally averaged approxi- 
mation by [oiiBa96l[l973), 



g ff (v, T) = 1.5 In (r/K) - In (v/GHz) - 5.307 



(5) 



For high frequencies, exact treatments show that gs — * 1 (e.g. 
Beckert et al. 2000) which we account for by a function that 
steps in when eq. (O would become < 1 . 

In Fig. |2] we show a free-free emission model fit to the ra- 
dio SED of NGC 1068. The model parameters (noted as model 
F) are shown in Table [3] We attenutated the model with ISM 
dust of optical depth Ty — 6 to be consistent with the NIR 
measurements, which concerns mainly the /-band data point 
in this case. The low-frequency cutoff occurs at the transition 
where the medium becomes optically thick to free-free emis- 
sion. With out parameters, we obtain a cutoff frequency of 
v(rff = 1) = 13.1 GHz. For the size of the emission region, we 
us ed the resolved core di ameter of 5 and 8.4 GHz observations 
bv lGallimore et all 12004, semi-major axis ~ 1 1 mas). These au- 
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thors also provide parameters for a free-free model, based on 
fluxes at their observed frequencies, which are approximately in 
agreement with our model using a broader frequency range. 



4. The influence of the radio components to the IR 

Since there is a gap between the 100 GHz and lOTHz regime, 
it is difficult to judge how the IR and radio components are 
connected. In particular, a high frequency turnover for the ra- 
dio component cannot be determined from the existing data. As 
a consequence, the radio source either breaks before the IR or 
reaches IR frequencies and, thus, may contribute to the unre- 
solved IR fluxes according to our models as discussed in the pre- 
vious section. The latter possibility leaves two open questions: 
(1) Is the unresolved flux in the long -baseline MIR interferomet- 
ric data from VLTI/MIDI coming from the central AGN, i.e. the 
radio source? (2) Can the radio source (especially a possible syn- 
chrotron source) be responsible for a signi ficant fraction of th e 
unresolved K band flux as observed by Wittkowski et al. (2004)? 
To address these questions, we extended the two investigated ra- 
dio models into the IR regime, as discussed in the previous sec- 
tion (Figs[TJ&[2] green dotted lines). 

We assume that the radio source is surrounded by the dusty 
torus. As a result, it will suffer from extinction by the dust within 
the torus (Figs[T]&[2] green dashed lines). The extiction that were 
used to reproduce the /-band observation with the radio emis- 
sion models are, however, much smaller than what is inferred 
from the lack o f optical and infrared broad lines {Ay ^ 50 mag 
iLutz et al]|2000h . Thus, the actual IR contribution of the radio 
emission source is most likely less than what is derived from our 
approach to reproduce the /-band data. This enforces the analy- 
sis in terms of obtaining strong upper limits for the contribution 
to the unresolved ^T-band and mid-IR flux. 

As can be seen from Figs [TJ & [2] the maximum possi- 
ble extrapolated mid-IR flux of the synchrotron and free-free 
power law is still a factor of 5-10 below the unresolved 78 m 
baseline MIDI flux. As a conservative limit, we estimate the 
maximum flux contribution of the radio source to the MIR 
FRadio(MIR)/F(MIR) <K 0.2 if it is sychrotron emission. This 
means that the unresolved MIDI flux is most likely originating 
exclusively from thermal dust emission in the torus rather than 
from the radio compon e nt. As shown in the previous section (see 
also iHonig et alj|2006l 12007). the unresolved emission can be 
fully interpreted in terms of a clumpy torus model. 

In the near-IR, the constraints are even tighter since (1) the 
assumption that 100% of the /-band flux is very conservative 
and (2) the difference in "true" (= broad-line region) and mod- 
elled extinction affects near-IR frequencies stronger than mid-IR 
frequencies. As can be seen in FigsQ]&[2] the non-torus contri- 
bution to the TT-band emission is certainly much less than 20%, 
no matter if the radio emission is due to synchrotron or free-free 
emission. This very strong upper limit is important for the inter- 
pretation of the unresolved flux measured by Wittkowski et al. 
(12004 using VLTI/VINCI: Even though the size of this un- 
resoved emissio n is < 0.3 pc, it must originate from the torus. 
Kishi moto el al] (120071) show that this size is consistent with 
the sublimation radius as expected from reverberation mapping 
of type 1 objects and assuming a central AGN luminosity of 
10 45 erg/s. As an alternative, substruct ure in a clu mpy torus can 
also account for the unresolved flux (Honig et al. 2007). 



5. Quantifying the contribution of the accretion disk 

An alternative interpretation for the unresolved VLTI/VINCI 
flux is contribution from the putative accretion disk in the center 
of NGC 1068, as discussed by IWitfkowski et al] (2004). In that 
case, the torus requires to be clumpy so that the accretion disk 
can be seen through "holes" of low or no extinction within the 
torus. 

In a clumpy medium with high optical depth of each dust 
cloud, the probability, Pad, that an accretion disk photon crosses 
all clouds without being absorbed can be estimated by Poisson 
statistics, Pa d - exp(-JV) (N — nu mber of clouds along the 
line-of-sight: iNatte & Panagial 119841) . From the edge-on torus 
model as presented in Sect. [3] we obtained N — 10. This re- 
sults in a probability to see the accretion disk through the torus 
clouds of Pad ~ 10~ 5 , which makes the scenario very unlikely. 
Nevertheless, we aim for constraining the possible accretion 
disk contribution to the A'-band flux. For that, we used a stan- 



dard accretion disk spectrum of the shape F Y 



,1/3 



exp(v/v max ) 



(e.g. JShakura & Sunvaevl[l973l) . The turnover frequency, v max , 
is characterised by the accretion disk temperature at the in- 
nermost stable orbit (~ 3PschwarzschiidX which can be inferred 
from the mass of the central black hole (Mbh ~ 10 7 Mq, 
Greenhill et al. 1996) and the modelled bolometric luminosity 
(L = 2 x 10 45 erg/s, see Sect.[3Tj- For NGC 1068, we obtain 
^max ~ 8 x 10 16 Hz which corresponds to a temperature of ap- 
proximately 1.5 x 10 6 K. We extended the accretion disk spec- 
trum down to IR wavelengths, so that the integrated accretion 
disk spectrum equals the bolometric luminosity as inferred from 
torus modelling. This allows to get an idea of the thermal emis- 
sion from the accretion disk, although it is completely obscured 
in the UV and optical. 

In Fig. [3] we show a comparison of the observed fluxes with 
the clumpy torus model (Sect. 13. It , and the accretion disk spec- 
trum which is extrapolated from the UV to the MIR (dotted line) 
and attenuated by dust extinction with Ty = 9.6 (dashed line) 
to reproduce the /-band flux. The obscured accretion disk spec- 
trum is obviously not able contribute significantly to the NIR or 
MIR emission. Again, it should be noted that the actual extinc- 
tion towards the broad line region is much higher than used to 
obtain the upper limit as described in the previous section. Thus, 
we safely conclude that that the unresolved /f-band flux as ob- 
served with VLTI/VINCI (40% of the total nuclear flux) does 
not originate from the accretion disk. 



6. Summary and Conclusions 

We analysed the SED of NGC 1068 from 10 9 Hz to 10 14 Hz. The 
SED was modelled by two components: (1) an IR component 
which is represented by a clumpy torus model, and (2) a ra- 
dio component represented either by synchrotron emission with 
F v oc v 1 / 3 , or optically thin free-free emission. For the infrared, 
we used our clumpy torus model to simultaneously reproduce 
the high-resolution NIR-to-MIR SED and the VLTI/MIDI cor- 
related fluxes at 78 m baseline. The model is characterized by a 
radial cloud density distribution oc r~ 3 ^ 2 and a half opening angle 
of 40°. The central AGN is obscured by 10 dust clouds on aver- 
age along a radial line of sight in the equatorial plane. For the ra- 
dio source, we showed that the power-law index of 0.29 + 0.07 is 
consistent with a spectral slope proportional to v 1 ^ 3 which is ex- 
pected from synchrotron emission by quasi-monoenergetic elec- 
trons. Further more, we showed that the radio emission can also 
originate from thermal free-free emission. Both radio models re- 
quire a low-frequency turnover between 1 and 10 GHz where 
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Fig. 3. Accretion Disk contribution to the IR spectrum. A standard accretion disk spectrum F v oc y 1 / 3 exp(v/v max ) was adjusted for 
the assumed parameters L — 2 x 10 45 erg/s and M BH ~ 10 7 Mq of NGC 1068 (green dotted line) and attenuated by dust obscuration 
(green dashed line) with ry = 9.6 to be consistent with the NIR measurements. 



the spectral slope changes. The steepness in the range between 
1-6 GHz suggests that the turnover is caused by a foreground 
screen of thermal plasma which becomes optically thick below 
~3-6GHz (F v oc exp(-(v/3 GHz) 21 )). This can be either con- 
nected to the gas in the galaxy or in the torus. 

For both scenarios - synchrotron and free-free emission - 
we presented models to extract physical parameters from the ra- 
dio SED, depending on the emission mechanism. The observed 
source sizes at 5 and 8.4 GHz are in agreement with synchrotron 
emission coming either from a disk (models A & B) or from 
small emission regions connected to the disk (e.g. magnetical 
reconnection zones in the disk corona) and possibly additional 
electron scattering in the surrounding medium (model C). The 
necessary Lorentz factors are y > 10 4 with B-field strengths of 
the order of 0.1 - 1 G. In case the radio SED is due to ther- 
mal free-free emission (model F), the observed fluxes and source 
sizes require an electron density n e = 8 x 10 5 cirT 3 and a tem- 
perature T = 1.3 x 10 6 K. 

The main goal of this study was to spectrally disentangle the 
near- and mid-infrared torus dust emission from possible con- 
tamination by the radio source. For that, we extrapolated the 
free-free and synchrotron models to the infrared. It was possible 
to show that the radio source is contributing certainly <20% to 
the observed unresolved flux at 78 m baseline. A similar anal- 
ysis shows that the unresolved ^T-band flux as observed with 
VLTI/VINCI is originating from the torus and almost uncontam- 
inated by the radio source. We were also able to argue against 
the possibility that in the 7^-band the accretion disk might be 
seen through holes in the clumpy torus. This means that all ob- 
served near- and mid-IR emission - resolved and unresolved - 



is thermal dust emission from the torus. This interpretation is 
consistent with the clumpy torus model that was presented. It 
is, however, not clear if the unresolved ^T-band flux represents 
emission from the inn er boundary of the torus as indicated by 
Kishimoto el al. (2007), or due to the clumpy substructure of the 
torus ( JHonigeT al. 2006; Honig & Beckeril2007l) . 
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